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Using primarily experimental inputs for S(Bd — > ipK s ), AMg s , AM_B d , BR (B — > tv) and tK 
along with necessary inputs from the lattice, we find that the measured value of sin(2/3) is smaller 
than expectations of the Standard Model by as much as 3.3 o, and also that the measured value 
of the BR(_B — > tv) seems to be less than the predicted value by about 2.8 a. However, through a 
critical study we show that most likely the dominant source of these deviations is in Bm s ) mixings 
and in sin(2/J) and less so in B — ¥ tv, and also that the bulk of the problem persists even if 
input from ei< is not used. The fact that kaon mixing and en are not the dominant source of the 
deviation from the Standard Model has the very important consequence that model independent 
considerations imply that the scale of the relevant new CP- violating physics is below 0(2 TeV), and 
possibly even a few hundred GeVs, thus suggesting that direct signals of the new particle(s) may 
well be accessible in collider experiments at the LHC and perhaps even at the Tevatron. 
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Introduction. It has been clear for quite sometime now 
that the CKM-paradigm p] of the Standard Model (SM) 
provides a quantitative description of the observed CP vi- 
olation, simultaneously in the B-system as well as in the 
K-system with a single CP-odd phase, to an accuracy of 
about 20% [2J. This major milestone in our understand- 
ing of CP violation was, of course, achieved in large part 
due to the spectacular performance of the two asymmet- 
ric B-factories. However, it should be recognized that in- 
put from the lattice for various weak matrix elements has 
also played a crucial role; in particular, the input from 
the lattice for the precise value of the hadronic parame- 
ter ("-Br-") characterizing the kaon-mixing amplitude is 
essential to demonstrate a quantitative understanding of 
the indirect CP-violation in Kt, — > tttt, i.e. ~ 10~ 3 . 
While the success of the CKM picture is very impressive, 
the flip side is that an accuracy of 0(20%) leaves open 
the possibility of quite sizable new physics contributions. 
Indeed, in the past few years as better data and better 
theoretical calculations became available some tensions 
have emerged [3- 8 J . It is clearly important to scruti- 
nize these tensions carefully to see if they persist and 
if so what they imply for the possible existence of new 
physics. 

In this paper we show that use of the latest experimen- 
tal inputs along with a careful use of the latest lattice 
results leads to a rather strong case in favor of a possible 
failure of the CKM picture for a sizable contribution due 
to beyond the Standard Model sources of CP violation 
that in sin 2/3 could be around 15-25%. Clearly if this 
result stands further scrutiny it would have widespread 
and significant repercussions for experiments at the 
intensity as well as the high energy frontier and, of 
course, also for our theoretical understanding. We will 
also show that we are able to isolate the presence of new 
physics to AB = 2 mixing amplitudes with a possible 
sub-dominant effect in kaon-mixing. In particular our 



analysis indicates that the data does not seem to provide 
a consistent interpretation for the presence of large new 
physics contribution to the tree amplitude for B — > tv. 

Choice of lattice inputs. Key inputs from experiment 
and from the lattice needed for our analysis are shown in 
Table [T] Below we briefly remark on a few of the items 
here that deserve special mention: 

• With regard to lattice inputs we want to emphasize that 
quantities used here have been carefully chosen and are 
extensively studied on the lattice for many years. Results 
that we use are obtained in full QCD with Np = 2 + 1 
simulations, are in the continuum limit, are fairly robust 
and emerge from the works of more than one collabora- 
tion and in most cases by many collaborations. 

• For using B — > tv in the context of our fits one clearly 
needs fg. Although this can be, and is directly calcu- 
lated on the lattice and results are known (see Table [I]), 
we will not use this as an input. Rather, our strategy 
is to extract the "predicted" or fitted value of Jb in the 
context of a given hypothesis so that comparison with the 
value directly determined from lattice calculations then 
can serve as a useful test of the viability of that partic- 
ular hypothesis. As we will show, this proves to be a 
rather useful criteria in checking the internal consistency 
of the hypotheses. This strategy will help us to isolate 
the main source of the deviation between experiment and 
the SM and thereby it may also serve as a useful way to 
uncover the nature of the underlying new physics. For 
completeness let us mention that fg needed for B — > tv 
in conjunction with our analysis is deduced indirectly by 
using the SU(3) breaking ratio, £ = f B ,Bl /2 / JbB 1 / 2 , 
f Bs Bl /2 and B d . 

• Regarding calculations of Bk on the lattice, it is useful 
to note that in the past 3 years a dramatic reduction in er- 
rors has been achieved and by now a number of indepen- 
dent calculations with errors < 5% and with consistent 
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Am Bd = (0.507 ± 0.005) ps" 1 a = (89.5 ± 4.3)° 

Am Bs = (17.77 ± 0.12) ps" 1 k s = 0.94 ± 0.02 [25H27] 

S^ks = 0.668 ± 0.023 [2SJ 7 =(74 ±11)° 

m c (m c ) = (1.268 ± 0.009) GeV B K = 0.740 ± 0.025 

m t , poie = (172.4 ± 1.2) GeV f K = (155.8 ± 1.7) MeV 

fB. V% = (276 ± 19) MeV [23] e K = (2.229 ± 0.012)10~ 3 

/s = (208 ± 8) MeV [23l [24] a B d = 1.26 ± 0.10 [H [24] 

Bb^tv = (1.68 ±0.31) x IP" 4 EHHSII 

a Our value of fg reflects the change in the overall scale (ri) 
recently adopted by the Fermilab/MILC and HPQCD collabora- 
tions ES] 

TABLE I: Inputs used in the fit. References to the original 
experimental and theoretical papers and the description of 
the averaging procedure can be found in Ref. [26]. Statistical 
and systematic errors are combined in quadrature. We adopt 
the averages of Ref. [26] for all quantities with the exception 
of \V u b\, £, fB s §l /2 , B K and f B (see text). 



central values have been obtained using Nf = 2 + 1 0- 
[T5] as well as Nf — 2 [2] dynamical simulations [T5] , 
Again, to be conservative we only use weighted average 
of two results that are both in full QCD, use different 
fermion discretizations (domain-wall and staggered) and 
that have no correlations between them [TO] [IT] . 

• Regarding V u b, since inclusive and exclusive methods 
differ appreciably, it should be clear that it is very diffi- 
cult to draw reliable conclusions by using this quantity; 
we will therefore make very limited and peripheral use of 
V u b only. For |V^b| e xci) to be conservative, we take the 
smaller of the two errors between the FNAL/MILC [T5] 
and HPQCD [T7] collaborations rather than taking the 
weighted average . Also, we add a 10% uncertainty to the 
inclusive determination of V u b in order to reflect intrin- 
sic uncertainties of the theoretical model adopted. The 
exclusive and inclusive determinations still differ at the 
1.8cr level; see Table [IJ 

• In a nutshell, we want to emphasize that we believe 
that the inputs used from the lattice are robust and 
reliable and therefore the implications resulting from 
their application should be taken seriously. 

Result of the Fit. The results of the fit are shown in 
Fig. [T] In Fig. [TJi we use as inputs from experiments, 
e K , AM s /AM d , 7 and BR(B -> tv) [33J and from the 
lattice, Bk, fg 3 B\/2 and Bd (but not fg) and we 



extract the fitted value of sin 2/3 and of fg. We obtain: 

sin(2/3) fit = 0.867 ± 0.048 , (1) 

which is about 3.3 a away from the experimentally mea- 
sured value of 0.668 ± 0.023. This is the key result of 
this paper providing a strong indication that the CKM 
description of the observed CP violation is breaking 
down [34] 

Clearly if this result is true it would be rather signif- 
icant; therefore we want to carefully check how robust 
it is and also try to isolate which of the relevant physi- 
cal process or processes may be seeing the effect of new 
physics. 

An important clue is provided by extracting the fitted 
value of fg along with the predicted value of sin(2/3) 
given above, we find: 

/|* = (199.7 ± 8.7) MeV . (2) 

This "predicted" value is in very good agreement with the 
one obtained by direct lattice calculation, fg = (208 ±8) 
MeV. This is a useful consistency check signifying that 
the SM description of the inputs used, especially of 
B — > tu, is working fairly well and that it is unlikely 
that the B — > tv tree amplitude is receiving large con- 
tributions from new physics. Most likely the dominant 
effect of new physics is in fact in sin(2/3). Note also that 
the fitted value of \V uh \ = (46.9 ± 3.7) x 10~ 4 is quite 
consistent with that determined by the inclusive method 
and deviates significantly (ss 3.6 a) from that obtained 
by the exclusive approach. Later we will comment on 
this regarding especially the issue of BR(B —> tv)). 

In order to further scrutinize the tentative conclusion 
reached above, we next present an alternate scenario de- 
picted in Fig.JlJi. Here, we make one important change in 
the inputs used. Instead of using the measured value of 
BR(_B — > tv) we now use as input the measured value of 
sin(2/3) from the "gold-plated" B^ — > i/jK s mode. Again, 
this fit yields two important predictions: 

BR(B -> rvf 1 = (0.754 ± 0.093) x 10~ 4 , (3) 
ff = (185.2 ± 8.6) MeV . (4) 

Eq. ([3]) deviates by 2.8er from the experimental mea- 
surement, as can also be gleaned from an inspection of 
Fig. [T|d. It is particularly interesting that also the fit 
prediction for fg now deviates by about 1.9tr from the 
direct lattice determination given in Table [TJ We believe 
this provides additional support that the measured value 
of sin(2/3) being used here as a key input is not consistent 
with the SM and in fact is receiving appreciable contri- 
butions from new physics. 

This leads us to conclude that while the presence of 
some sub-dominant contribution of new physics in B — > 
tv is possible, a large contribution of new physics in there 
is not able to explain, in a consistent fashion, the tension 
we are observing in the unitarity triangle fit. 
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FIG. 1: Unitarity triangle fit. In the top panel, the contour 
and the fit predictions for sin(2/3), fs and \V u b\ are obtained 
using Vcb, sk, 7, AMs d s and B — s> In the bottom panel, 
the contour and the fit predictions for BR(B — > tv), /b and 
\V u b\ are obtained using V&, £k, 7, AMs d s and sin(2/3). 



This conclusion is corroborated by the observation that 
even without using B — > tv at all, and using as input only 
e K , AM B JAM Bd and |F c6 | (see Fig. [I]), the predicted 
value of sin(2/3) deviates by 2.1a from its measurement 
(in this case we find sin(2/3) fit = 0.829 ± 0.079). Thus, 
possible new physics in B — > tv can alleviate but not 
remove completely the tension in the fit. 

We recall that the fit above is actually the simple fit we 
had reported some time ago (now with updated lattice 
inputs) with its resulting « 2 it deviation [1]. This fit is 
somewhat special as primarily one is only using AF = 2 
box graphs from ex and AMbJ AMs d in conjunction 
with lattice inputs for Br and the SU(3) breaking ratio 
£. The experimental input from box graphs is clearly 
short-distance dominated and for the lattice these 
two inputs are particularly simple to calculate as the 
relevant 4-quark operators have no mixing with lower 
dimensional operators and also require no momentum 
injection. The prospects for further improvements in 
these calculations are high and the method should 
continue to provide an accurate and clean "prediction" 
for sin(2/3) in the SM. So even if the current tensions 
get resolved, this type of fit should remain a viable way 
to test the SM as lattice calculations and experimental 
inputs continue to improve. 
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FIG. 2: Unitarity triangle fit without semileptonic decays (up- 
per panel) and without use of K mixing (lower panel). 



Roles of V c i,, Ek, V u b and of hadronic uncertain- 
ties. The fit described above does use V c h where again 
the inclusive and exclusive methods differ mildly (about 
1.7a). Of greater concern here is that ex scales as \V c b\ A 
and therefore is very sensitive to the error on V c b- We 
address this in two ways. First in Fig. [2^ we study a fit 
wherein no semi-leptonic input from b — > c or b — > u is 
being used. Instead, in this fit BR(v3 — > tv) and AMb b 
along with ek, AMbJ AMB d and 7 are used. Interest- 
ingly this fit gives 

sin(2/3) fit = 0.891 ± 0.052 , (5) 
/|* = (200.7 ± 8.6) MeV (6) 

Thus, once again, sin(2/3) is off by 2.8<r whereas Jb is in 
very good agreement with directly measured value which 
we again take to mean that the bulk of the discrepancy 
is in sin(2/3) rather than in B — > tv or in V c b. 

Next we investigate the role of ex- In Fig. [2Jd we 
show a fit where only input from B-physics, namely 
AM B JAM Bd , AMb b , 7, \Vcb\ and BR(S -> tv) are 
used. This fit yields, 

sin(2/3) fit = 0.891 ± 0.054 , (7) 
/|* = (195 ±11) MeV (8) 

Thus, sin(2/3) fit is off by w 2.4er and again /|* is in good 
agreement with its direct determination. We are, there- 
fore, led to conclude that the role of ck in the discrepancy 
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FIG. 3: Unitarity triangle fit with V u b- See the caption in 
Fig-0 



is subdominant and that the bulk of the new physics con- 
tribution is likely to be in B-physics. As before, the fact 
that the fitted value of fg is in good agreement with 
its direct determination seems to suggest that the input 
BR(/3 — > tv) is most likely not in any large conflict with 
the SM, though, obviously we cannot rule out the pos- 
sibility of it receiving a sub-dominant contribution from 
new physics. 

Although we believe that we have been very careful 
in taking the input from lattice calculations and their 
associated errors (see Table to gain further confidence 
we study the effect of increasing the total error in each 
of the input quantity by 50%, we find that qualitatively 
little change takes place from Eq. ([I]): 

sin(2/3) fit = 0.854 ± 0.052 , (9) 
fs = (202 ± 13) MeV (10) 

again we find a 3.0 a deviation in sin(2/3) from the mea- 
sured value and the fitted value of f& in very good agree- 
ment with the direct lattice determination strongly sug- 
gesting once again that the discrepancy with the CKM 
is rather serious. 

So far we have stayed away from using \V u b\ because 
of the large associated uncertainties with it. As another 
rough consistency check, let us mention that with the in- 
clusion of \V u f,\ (resulting from combining inclusive and 
exclusive methods with an estimated error of about 15% 
- see Table IB, the results outlined above do not change 



FIG. 4: Summary of sin(2/3) determinations. The entry 
marked *** (ninth from the top) is obtained with lattice er- 
rors increased by 50% over those given in Table [I] for each of 
the input quantities that we use and the entry marked +++ 
(tenth from the top) corresponds to adding an hadronic un- 
certainty SAS^k = 0.021 to the relation between sin(2/3) and 
S^k- See the text for further explanations. 



qualitatively. The discrepancy between the fitted and 
measured values of sin(2/3) is mildly reduced to about 
2.5<t (see Fig. [3^i) while the fit result for fg is quite com- 
patible with its direct lattice determination. In Fig. [SJd 
we investigate the alternate hypothesis of using S(ipKs) 
as a SM input and find that even with the use of \V u b\ the 
results remain essentially unchanged: both BR(B — > tv) 
and fs deviate appreciably from their respective direct 
determinations . 

A compilation of all ten fits that we study for sin 2/3 
are shown in Fig. [I] Notice that there is only one case in 
here (8th from the top) where the discrepancy in sin 2/3 
is only 0(la). We believe this is primarily a reflection of 
the large (as 14.4%) uncertainty with our combined V u b 
fit originating from the large disparity between inclusive 
and exclusive determinations; this is why we make only 
a limited use of V u b in our fits. 

Now with regard to B — > tv, Fig. [5] shows a sum- 
mary of predictions versus the measured BR. Notice that 
whenever the measured value of sin(2/3) is used as an in- 
put, the predicted BR is » 2.8cr from the measured one. 
In the preceding discussion we have emphasized that this 
seems to us to be a consequence of new physics largely in 
B mixings. This conclusion receives further strong sup- 
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port when we try determine the B — tv branching ratio 
without using sin 2/3; indeed as shown in Fig. [5] when 
we use 6k, AMb 9 , V c b and 7 only the fitted value of 
BR(B — > tv) is in very good agreement with the mea- 
sured value. 

In principle, of course the prediction for BR(S — > tv) 
only needs the values of fg and of V u f,. Fixing now 
fs = 208 ± 8 MeV as directly determined on the lat- 
tice (see Table [I]) we show the corresponding two pre- 
dictions for the BR using separately the values of V u b 
determined in inclusive and exclusive decays. It is clear 
that the inclusive determination yields results that are 
within one a of experiment (see also Fig. [I]); however 
with V u b from exclusive modes (that makes use of the 
semileptonic form factor as determined on the lattice), 
the BR deviates by ~ 3<r from experiment. This may 
be a hint that lattice based exclusive methods have some 
intrinsic difficulty or that the exclusive modes are sen- 
sitive to some new physics that the inclusive modes are 
insensitive to, e.g. right-handed currents [351 El]- In ei- 
ther case, this reasoning suggests that we try using the 
value of V u b given by inclusive methods only in our fit 
for determining sin 2/3. Wc then obtain results that are 
very compatible with the no-V u b case: the fitted value 
of sin(2/3) deviates by 3.2a from its direct determination 
(see Fig Ql . 

Finally, we address the possible presence of sizable 
hadronic uncertainties in the relation between sin(2/3) 
and Srf,K- Naive estimates of the impact that the CKM 
suppressed u-penguin amplitude, which is causing the 
"pollution", has on S(B —> J/tpK) point to a sub- 



percent effect. Quantitative studies based on QCD- 
factorization [37], perturbative QCD [38] and conserva- 
tive model independent bounds on possible rescatter- 
ing effects [35] corroborates [ID] the above mentioned 
naive expectations. An alternative approach based on 
use of flavor SU(3) to connect B — » J/tpK and B — » 
J/ipTt hadronic matrix elements has been proposed in 
Refs. [41] |42] . This method is based on the observation 
that, up to SU(3) corrections, the tree and penguin ma- 
trix elements appearing in B — > J /%pK and B — ¥ J/ipTr 
decays are identical. Since the penguin topology in 
B J/ipn is not CKM suppressed with the respect to 
the corresponding tree amplitude, time-dependent CP 
asymmetries in B — > J/i/jtt are highly sensitive to effects 
that affect the J/ipK mode at the percent level. Un- 
fortunately, data on B — > J/ipir CP asymmetries is not 
precise enough to offer a measurement of the penguin to 
tree ratio and phase and we do not reliably know how 
large should be the SU(3) breaking. The upper limits 
are presently more than an order of magnitude above 
all the direct estimates. For these reasons we believe 
that, presently, it is not useful to adopt B — > J/tpir de- 
cays as a sole handle on the size of penguin pollution 
in B — > J/tpK. For completion we mention that even 
adopting the estimate of Ref. [43] for penguin pollution 
in the extraction of sin(2/3), i.e. 5AS = 0.021, Eq. Q 
deviates from the measurement at the 3.2tr level. 
Summary, perspective & outlook. In this paper we 
have mainly emphasized that our analysis strongly sug- 
gests that the SM predicted value of sin(2/3) is around 
0.85 whereas the value measured experimentally via the 
gold plated tJjK s mode is around 0.66 constituting a de- 
viation of about 3 a from the SM (see Fig. |4| . To put 
this result in a broader perspective let us now recall that 
in fact in the SM sin(2/3) can also be measured via the 
penguin dominated modes (see Fig. [4]) [44H47] . Unfortu- 
nately these modes suffer from a potentially large pen- 
guin pollution, though there are good reasons to believe 
that the i]'K s , <\>K 8 and 3 K s modes are rather clean [4"8T - 
150] wherein the deviations from sin 2/3 are expected to be 
only 0(few %). The striking aspect of these three clean 
modes as well as many others penguin dominated modes 
(see Fig. [2]) is that the central values of almost all of them 
tends to be even smaller than the value (0.66), measured 
in ipK s , and consequently tend to exhibit even a larger 
deviation from the SM prediction of around 0.85. Thus, 
seen in the light of our analysis, the deviation in these 
penguin modes suggests the presence of new CP- violating 
physics not just in B- mixing but also in b — > s penguin 
transitions. 

Moreover, the large difference (ra (14.4 ± 2.9)%) [5T] 
in the direct CP asymmetry measured in B — > K + ir~ 
versus that in B + ^ K + ir° provides another hint that 
b — > s penguin transitions may be receiving the contri- 
bution from a beyond the SM source of CP-violation. 
To briefly recapitulate, in the SM one naively expects 
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this difference to be vanishingly small and careful es- 
timates based on QCD factorization ideas suggest that 
it is very difficult to get a difference much larger than 
(2.2 ±2.4)% [5J. 

Of course, if b — > s penguin transitions (AFlavor = 1) 
are receiving contributions from new physics, then it is 
quite unnatural for B s mixing amplitudes ( AFlavor = 2) 
to remain unaffected. Therefore, this reasoning suggests 
that we should expect non-vanishing CP asymmetries in 
B s — > ifxp as well as a non- vanishing di-lepton asymme- 
try in B s — > X s lv. As is well known, at Fermilab, in 
the past couple of years CDF and DO experiments have 
been studying CP asymmetry in B s — > ip(j>. The lat- 
est analysis with about 6.1 fb _1 luminosity at DO shows 
a deviation from the SM of about 2 a [52], which is a 
mild increase from their previous result (~ 1.8tr) with 
2.8 fb^ 1 whereas at CDF the results with corresponding 
increase in luminosity has shown a downward shift in the 
deviation from the SM to ps 0.8ct from the earlier s» 1.8a 
deviation |53j . Thus, the Fermilab findings in B s — > ip(j> 
are not yet clear. 

Another interesting and potentially very important de- 
velopment with regard to non-standard CP in B s is that 
in the last several months DO has announced the observa- 
tion of a large dimuon asymmetry in B-decays amounting 
to a deviation of (« 3.2ct) from the miniscule asymmetry 
predicted in the SM [5H [55J . They attribute this largely 
to originate from B s mixing. While this is a very excit- 
ing development, their experimental analysis is extremely 
challenging and a confirmation is highly desirable before 
their findings can be safely assumed. 

Be that as it may, we reiterate that our analysis sug- 
gests that the deviation from the SM in sin(2/3) is diffi- 
cult to reconcile with errors in the inputs from the lat- 
tice that we use, and strongly suggests the presence of 
a non-standard source of CP violation largely in B/B s 
mixings, thereby predicting that non-standard signals of 
CP violation in S(Bd — > rj'K s , 4>K S , 3K s etc.) as well as in 
S(B S — > ip4>), and the semileptonic and di-lepton asym- 
metries in B s , and possibly also in Bd, decays will persist 
and survive further scrutiny in experiments at the inten- 
sity frontiers such as Fermilab (CDF, DO), LHCb and the 
Super-B factories. Lastly, the fact that our analysis rules 
out the presence of large new physics in kaon mixing has 
very important repercussions for the mass scale of the 
underlying new physics contributing to these anomalies 
in B, B s decays: model independent analysis then imply 
that the relevant mass scale of the new physics is nec- 
essarily relatively low, i.e. below 0(2 TeV) and may, in 
principle, be as low as a few hundred GeVs [5J[5S]. Thus, 
collider experiments at the high energy frontier at LHC 
and possibly even at Fermilab should see direct signals 
of the underlying degrees of freedom appearing in any 
relevant beyond the Standard Model scenario. 
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